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Abstract 

Differential dijet cross sections have been measured with the ZEUS detector for pho- 
toproduction events in which the hadronic final state containing the jets is separated with 
respect to the outgoing proton direction by a large rapidity gap. The cross section has 
been measured as a function of the fraction of the photon (x® BS ) and pomeron {f3 OBS ) 
momentum participating in the production of the dijet system. The observed x® BS depen- 
dence shows evidence for the presence of a resolved- as well as a direct-photon component. 
The measured cross section da/df3 OBS increases as [3 OBS increases indicating that there 
is a sizeable contribution to dijet production from those events in which a large fraction of 
the pomeron momentum participates in the hard scattering. These cross sections and the 
ZEUS measurements of the diffractive structure function can be described by calculations 
based on parton densities in the pomeron which evolve according to the QCD evolution 
equations and include a substantial hard momentum component of gluons in the pomeron. 
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1 Introduction 



A successful phenomenological description of the available data on soft diffractive processes has 
been obtained using Regge theory and the exchange of a trajectory with the quantum numbers 
of the vacuum, the pomeron (IP). Hard processes in diffractive reactions |], HEIIl, H>EI0>0 
provide a tool to investigate the partonic nature of this colour- singlet exchange and to test the 
universality of its properties in different reactions. 

The first experimental evidence pointing to a partonic nature of the pomeron was the ob- 
servation of jet production in pp collisions with a tagged leading proton (or antiproton) made 
by the UA8 Collaboration || following the proposal of Ingelman and Schlein [0]. At HERA, 
the ZEUS and HI Collaborations observed neutral-current deep-inelastic ep scattering (DIS) 
events at high Q 2 (where Q 2 is the virtuality of the exchanged photon) characterised by a large 
rapidity gap from the proton direction |T^,|rTJ. The properties of these events were suggestive 
of a diffractive interaction mediated by pomeron exchange between a highly-virtual photon and 
a proton. The observed Q 2 dependence indicated a pointlike nature of the interaction and a 
leading-twist mechanism. Measurements of the diffractive structure function in DIS |~2[[l3|.|n 



were found to be consistent within the experimental uncertainties with a diffractive structure 
function which factorises into a pomeron flux factor, which depends on the momentum frac- 
tion lost by the proton (x/p), and a pomeron structure function, which depends on Q 2 and 
the momentum fraction of the struck quark within the pomeron (f3). The pomeron structure 



function showed an approximate scaling with Q at fixed (3. More recent HI measurements [15 



have been analysed in terms of parton densities in the pomeron which evolve according to the 
DGLAP equations |I6 . The observed scaling violations indicate that most of the momentum 



of the pomeron is carried by gluons. 

Diffractive photoproduction (Q 2 ph 0) of high transverse energy jets at HERA provides a 
process which is sensitive both to the quark (e.g. via — > qg) and gluon (e.g. via jg — > 
qq) densities in the pomeron. The results on the diffractive structure function in DIS (13 



combined with the measured inclusive jet cross sections in diffractive photoproduction gave 
the first experimental evidence for a gluon content of the pomeron [I7[|. The data indicated 
that between 30% and 80% of the momentum of the pomeron carried by partons is due to 
hard gluons. The conclusion was based on the assumption that the same pomeron flux factor 
( "Regge factorisation" ) and the same pomeron parton densities apply to diffractive DIS and to 
diffractive jet photoproduction at similar hard scales ("hard-scattering factorisation"). 



The CDF Collaboration has recently observed diffractive production of both iy-bosons [IS 



and dijet events [|I^| in pp collisions at ^fs = 1.8 TeV. The analysis of these measurements 
yields an estimate of the hard-gluon content of the pomeron of (70 ± 20)% [pl|, which is in 
agreement with the result obtained at HERA [|15,|l7|]. However, the fraction of the total pomeron 
momentum carried by partons, assuming the pomeron flux factor of Donnachie and Landshoff 



is well below our result [17]. This discrepancy could be an indication of a breakdown of 



hard-scattering factorisation for diffractive hard processes in pp interactions ||T9],pQ,|21| . 

Diffractive photoproduction of dijets is sensitive to the underlying two-body processes. In 
leading-order (LO) QCD two types of processes contribute to jet production |22|,|23[: either the 
photon interacts directly with a parton in the pomeron (the direct process) or the photon acts 
as a source of partons which scatter off those in the pomeron (the resolved process) . Examples 
of Feynman diagrams for diffractive dijet photoproduction are shown in Figure [|. The cross 
section dependence on the fraction of the photon momentum participating in the production 
of the two jets, x® BS , is sensitive to the presence of these two components p4|,[2 



1 



In this paper, new measurements of differential cross sections for diffractive dijet photopro- 
duction at HERA are presented. The results are presented as a function of the pseudorapidityQ 
[rf et ^j and transverse energy of each jet (E^*), the 7p centre-of-mass energy (W), x° BS and the 
fraction of the pomeron momentum (j3 OBS ) participating in the production of the two jets with 
highest Ej? in an event. These measurements together with those of the diffractive structure 
function [13] are analysed in terms of parton densities in the pomeron which are assumed to 
evolve according to the DGLAP equations. We therefore test whether it is possible to describe 
both sets of data with the same pomeron parton densities, assuming Regge factorisation. The 
data sample used in this analysis was collected with the ZEUS detector in e + p interactions at 
the HERA collider and corresponds to an integrated luminosity of 2.65 pb _1 . 



2 Kinematics of diffractive hard processes 

Diffractive photon-dissociative processes^ in e + p collisions are characterised by a final state 
consisting of a hadronic system X, the scattered positron and the scattered proton 

e+(k) + Pi {P) -> e + (k') + X + pf(P'), 

where pt (p/) denotes the initial (final) state proton. The kinematics of this process are described 
in terms of four variables. Two of them describe the positron-photon vertex and can be taken 
to be the virtuality of the exchanged photon (Q 2 ) and the inelasticity variable y defined by 

, 2 , .. /'•'/ 



Q 2 = -q 2 = -(k-k'Y and y 



P-k' 



The other two variables describe the proton vertex: the fraction of the momentum of the initial 
proton carried by the pomeron (xjp), and the square of the momentum transfer (t) between the 
initial and final state proton, defined by 

Xp = (P ~ p/) " 9 and t=(P- P'f. 
P ■ q 

The hard scale is given by Q in diffractive DIS and by the jet transverse energy in diffractive 
photoproduction (Q 2 « 0). In DIS, an additional and useful kinematic variable is given by 

,3 Q2 Q1 



2(P-P')-q syxp 

where y/s is the e + p centre-of-mass energy. In models where the pomeron has a partonic 
structure, (3 is the equivalent of the Bjorken-a; variable in e + lP interactions. 

In diffractive dijet photoproduction processes the hadronic system X contains at least two 
jets, 

e + (k) + p t (P)^ e + (k') + X + p f (P') -> e +(A/) + X(jet+jet + X r ) + p f (P'), (1) 



lr rhe ZEUS coordinate system is defined as right-handed with the Z axis pointing in the proton beam 
direction, hereafter referred to as forward, and the X axis horizontal, pointing towards the centre of HERA. 
The pseudorapidity is defined as rj = — ln(tan | ) , where the polar angle 9 is taken with respect to the proton 
beam direction. 

2 These processes will be henceforth referred to as diffractive processes unless otherwise stated. 



2 



where X r denotes the remainder of the final state X not assigned to the two jets with the 
highest E^ 1 . In two-to-two massless-parton scattering, the fractions of the photon (x 7 ) and 
pomeron (/3p) momenta carried by the initial-state partons (in the infinite-momentum frame 
of the parent particles) are defined by 

_ {pji + Pji) • q' , o _ {PJi+PJ2)-q 

Xry ailU , 

q ■ q' q ■ q 

where pji and pj2 are the momenta of the two final state partons, q 1 is defined as q' = P — P' 
and the approximations q 2 ~ q' 2 ~ have been used. For LO direct processes, x. 
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Since the partons are not measurable objects, the observables x° BS [25j and f3 OBS , defined 
in terms of jets, are introduced, 



x OBS _ l^jets-^T e an( ^ /JOBS _ l^jets -^T e 

7 2y£ e ^ 2x F E p ' 

where £? e is the incident positron (proton) energy and the sum runs over the two jets of 
highest E^ 1 in an event. The variable x® BS (f3 OBS ) is an estimator of the fraction of the photon 
(pomeron) momentum participating in the production of the two jets with highest E^. The 
LO direct and resolved processes populate different regions of x° BS , with the direct processes 
being concentrated at high values of x° BS . 

Diffractive processes give rise to a large rapidity gap between the hadronic system X and 
the scattered proton: Ay G ap = y Pf — Ymaxi where y Pf is the rapidity of the scattered proton 
and is the rapidity of the most-forward-going hadron belonging to the system X. Instead 
of y" x the pseudorapidity (^at) °f the most-forward-going hadron in the detector was used 
to select diffractive events. The same signature is expected for double dissociation where the 
scattered proton is replaced by a low-mass baryonic system (N). In this measurement, the 
outgoing proton (or system N) was not observed. 



3 A model for diffractive hard processes 

Several models have been developed to describe diffractive hard processes assuming that the 
main mechanism for reaction ([!]) is pomeron emission by the proton. The measurements of 
the diffractive structure function [[L3|] and of the dijet cross sections in diffractive photoproduc- 



tion presented here are analysed in terms of a model in which both Regge and hard-scattering 
factorisation are assumed ^, [| [|, || (factorisable model). In this model, first proposed by 
Ingelman and Schlein resolved- and direct-photon contributions to dijet diffractive pho- 
toproduction can be calculated. Additional coherent-pomeron contributions are expected in 
resolved photoproduction from processes in which the whole pomeron initiates the hard scat- 
tering || 0,0. In principle such contributions would lead to factorisation breaking. However, 
the inclusion of a delta-function term at (5 — 1 in the gluon density of the pomeron in the 
factorisable model would lead to a similar (3 distribution. 

In QCD, the factorisation theorem ensures that the parton densities (for a given hadron) 
extracted from measurements of DIS apply also to other inclusive hard processes involving 
the same initial-state hadron. This theorem has been recently extended to diffractive DIS 



27 which justifies our analysis of the leading-twist diffractive structure function in terms of 



parton densities in the pomeron. This new theorem does not address the question of Regge 
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factorisation, which relates the properties of the pomeron in diffractive DIS to those measured 
in hadron-hadron reactions. 

The proof of hard-scattering factorisation fails for diffractive hard processes in hadron-hadron 
collisions ||27|| . Therefore, the contribution of the resolved process to diffractive dijet photopro- 



duction is expected to be non-factorisable. In what follows, hard-scattering factorisation is 
assumed to hold for both processes (direct and resolved) in the range x 7 > 0.2. We therefore 
check whether a consistent set of pomeron parton densities is able to describe simultaneously the 
measurements of the diffractive structure function in DIS and the cross sections for diffractive 
dijet photoproduction within the precision of the present data. 

To summarise, in this factorisable model, the pomeron is assumed to be a source of partons 
which interact either with the photon (in DIS and in the direct process of photoproduction) 
or with a partonic constituent of the photon (in the resolved process of photoproduction). 
Calculations based on such a model involve three basic ingredients: a) the flux of pomerons 
from the proton as a function of xp and t, b) the parton densities in the pomeron, and c) 
the matrix elements for the hard subprocess. For the resolved process in photoproduction, a 
fourth ingredient is needed: the parton densities in the photon. The pomeron flux factor is 
extracted from hadron-hadron collisions using Regge phenomenology, and the matrix elements 
are computed in perturbative QCD. However, the parton densities are a priori unknown and 
have to be extracted from experiment. 

In what follows, the pomeron flux factor given by Donnachie and Landshoff (DL) [[| has 
been used, 

/ I +\ ^0 2 l_2a(t) 

fp/p{X]p, t) = — Fi (t) xjp , 

where Fi(t) is the elastic form factor of the proton, 6 the pomeron-quark coupling (6 = 
1.8 GeV" 1 ) and a(t) the pomeron trajectory (a(t) = 1.085 + 0.25t with t in GeV 2 ) taken from 
hadron-hadron data. 

The parton densities in the pomeron, fi/jp({3, fi 2 ) with i = q,g, depend upon the fraction (3 
of the pomeron momentum carried by parton i and the scale fi at which the parton structure 
of the pomeron is probed. In diffractive dijet photoproduction processes, one choice, which 
is used here, is /i = px, where px is the transverse momentum of either of the two outgoing 
partons. The scale [i is set equal to Q for diffractive DIS. The parton densities in the pomeron 
are evolved in /j, 2 according to the DGLAP equations. 

In this analysis (see Section 9), the parton densities in the pomeron are determined from 
a simultaneous fit to the ZEUS measurements of the diffractive structure function F 2 D (/3,Q 2 ) 
(see below) |13|] and the measurements of diffractive dijet cross sections in photoproduction 
presented in Section 8. The shapes of the resulting parton densities do not depend on the 
normalisation of the pomeron flux factor and depend only weakly on the xjp-functional form 
as long as this is the same in diffractive DIS and dijet photoproduction. 

3.1 The diffractive structure function 

For unpolarised beams, the differential cross section for diffractive DIS can be described in 
terms of the diffractive structure function F®^\[3,Q 2 ,X]p,t) 



dpdQ 2 dxipdt j3Q A 

where a is the electromagnetic coupling constant and the contribution of the longitudinal 
diffractive structure function is neglected. An integration over the entire range of t defines the 
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diffractive structure function F 2 D ^ 3 \(3,Q 2 ,X]p) as measured in [13 



~ — (l + (1 - y?)F^\^Q\ X]P ). 



d(3dQ 2 dxip (3Q A 

To illustrate the (3 and Q 2 dependence of F 2 D< ^\(3, Q 2 , x P ), the structure function F 2 D (f3,Q 2 
was also measured in 



F 2 D (f3,Q 2 ) = r^dxjpF^iP^.xjp) 



In the factorisable model, the diffractive structure function F 2 ^ is assumed to be of the 
form 



F 2 D{A \p,Q 2 , X]P ,t) = f P/p (x P ,t)-FF(P,Q 2 ), 



where at LO the pomeron structure function F 2 depends on the pomeron parton densities as 
given by 

F 2 F ((3,Q 2 ) = Y,epU/AP,Q 2 )- 

j 

In this expression, ej is the electric charge of quark qj and the sum runs over all quark flavours 
which contribute at the given value of Q 2 . Therefore, the LO calculation for F 2 in this model 
is of the form 



F 2 D (P,Q 2 ) =Y, e Pf*MP>Q 2 ) ■ r^dxjp dt fjp /p {xjp,t). 

•'ZjPmin 



At next-to-leading order (NLO), F 2 also depends on the gluon density in the pomeron. The 
measurements of the diffractive structure function are analysed (see Section 9) in terms of NLO 
QCD calculations. 



3.2 Dijet cross sections in diffractive photoproduction 

The dijet cross sections in diffractive photoproduction contain contributions from both the 
direct and resolved processes. As an example, the contribution of the direct process to the 
cross section for reaction (|]) is given by 

Vdir = J dyf j/e (y) J J dx P dtft>/ p (x P ,t)^2 J d/3j2 J d Pr d<Tt+ ^ 3+k (g, jj 2 ) ■ f i/P (/3, /j 2 ), 



where / 7 / e is the flux of photons from the positronP]. The sum in i runs over all possible types 
of partons present in the pomeron. The sum in j and k runs over all possible types of final state 
partons and ai +1 ^j + k is the cross section for the two-body collision i + 7 — >• j + k and depends 
on the square of the centre-of-mass energy (s), the transverse momentum of the two outgoing 
partons (pr) and the momentum scale at which the strong coupling constant (a s (/j, 2 )) is 
evaluated. 



The Q 2 dependence has been integrated out using the Weizsacker- Williams approximation. 
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4 Experimental conditions 



During 1994 HERA operated with protons of energy E p = 820 GeV and positrons of energy 
E e = 27.5 GeV. The ZEUS detector is described in detail in |29],|30|]. The main subdetectors 
used in the present analysis are the central tracking system positioned in a 1.43 T solenoidal 
magnetic field and the uranium-scintillator sampling calorimeter (CAL). The tracking system 
was used to establish an interaction vertex and to cross-check the energy scale of the CAL. 
Energy deposits in the CAL were used to find jets and to measure jet energies. The CAL is 
hermetic and consists of 5918 cells each read out by two photomultiplier tubes. Under test beam 
conditions, the CAL has energy resolutions of 18%/\/E for electrons and 35%/y/E for hadrons. 
Jet energies are corrected for the energy lost in inactive material in front of the CAL. This 
material is typically about one radiation length. The effects of uranium noise were minimised 
by discarding cells in the inner (electromagnetic) or outer (hadronic) sections if they had energy 
deposits of less than 60 MeV or 110 MeV, respectively. The luminosity was measured from 
the rate of the bremsstrahlung process e + p — > e + p~y. A three-level trigger was used to select 
events online [Q. At the third level, where the full event information is available, the events 



were required to have at least two jets with jet transverse energy in excess of 3.5 GeV and jet 
pseudorapidity below 2.0 reconstructed using the CAL cell energies and positions as input to a 
cone algorithm (see Section 6). 



5 Monte Carlo simulation 

The response of the detector to jets and the correction factors for the cross sections for dijet 
production with a large rapidity gap were determined from Monte Carlo samples of events. 

The program PYTHIA 5.7 [|3l[] was used to generate standard (non-diffractive) hard pho- 
toproduction events for resolved and direct processes. The photon momentum spectrum was 
calculated using the Weizsacker- Williams approximation. Events were generated using GRV- 
HO |32| for the photon parton distributions and MRSA [33] for the proton parton distributions. 



The partonic processes were simulated using LO matrix elements, with the inclusion of initial- 
and final-state parton showers. Fragmentation into hadrons was performed using the LUND 



string model [|34| as implemented in JETSET ||35|| . Samples of events were generated with 
different values of the cutoff on the transverse momentum of the two outgoing partons, starting 
at pTmin = 2.5 GeV. 



Diffractive processes were simulated using the program POMPYT 2.5[] [35|. This is a Monte 
Carlo program where, within the framework provided by PYTHIA, the proton emits a pomeron 
whose partonic constituents subsequently take part in a hard scattering process with the pho- 
ton or its constituents. For the resolved processes, the parton densities of the photon were 
parametrised according to GS-HO |37|]f] evaluated at p T . The parton densities in the pomeron 
were parametrised according to a hard distribution (3(1-/3) and the DL form was used for the 
pomeron flux factor. 

All generated events were passed through the ZEUS detector and trigger simulation programs 



30 1 . They were reconstructed and analysed by the same program chain as the data. The 



resulting Monte Carlo distributions agree reasonably well with the data. 



4 This version of POMPYT has been modified to make use of pomeron parton densities which evolve with 
the scale according to the DGLAP equations. 

5 The correction factors depend very weakly on the specific set of photon parton distributions used. 
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6 Jet search and reconstruction of kinematic variables 



An iterative cone algorithm in the n — <p plane is used to reconstruct jets from the energy 
measured in the CAL cells for both data and simulated events, and also from the final-state 
hadrons for simulated events. A detailed description of the algorithm can be found in [ 38[ . The 



jets reconstructed from the CAL cell energies are called cal jets and the variables associated 
with them are denoted by E^ l cal , rf c ^ and ip 3 ^. The axis of the jet is defined according to the 

Snowmass convention |35| , where n 3 ^ ((fifi) is the transverse-energy weighted mean pseudora- 



pidity (azimuth) of all the CAL cells belonging to that jet. The cone radius used in the jet 
search was set equal to 1. 

For the Monte Carlo events, the same jet algorithm is also applied to the final-state particles. 
The jets found are called hadron jets and the variables associated with them are denoted by 
E T%d, rfhad, and Vhad- Hadron jets with E^ had > 6 GeV and -1.5 < rf h G * d < 1 are selected. 

The comparison of the reconstructed jet variables between the hadron and the cal jets in 
simulated events |^0[ shows no significant systematic shift in the angular variables rf c ^ and 
Pcai w hh respect to nj^ d and pj^d- The resolutions are 0.07 units in n 3 ,^ and 5° in (p 3 ^. The 
transverse energy of the cal jet underestimates that of the hadron jet by an average amount 
of 16% with an r.m.s. of 11%. The transverse energy corrections to cal jets averaged over the 
azimuthal angle were determined using the Monte Carlo samples of events. These corrections 
are constructed as multiplicative factors, C(E 3 rl al} n 3 ca i l ), which, when applied to the Et of the 
cal jets, give the true transverse energies of the jets, E^ = CiE^^rf^) x E^l al These 



corrections mainly take into account the energy losses due to the inactive material in front of 
the CAL. 



The 7p centre-of-mass energy W = ^Jys is estimated using the method of Jacquet-Blondel HT 



Wj B = y/2E p (E-P z ), (2) 

where E is the total energy as measured by the CAL, E = J2iEi, and Pz is the Z-component 
of the vector P = Y^iE^; in both cases the sum runs over all CAL cells, Ei is the energy of 
the calorimeter cell i and r*j is a unit vector along the line joining the reconstructed vertex and 
the geometric centre of the cell i. Due to the energy lost in the inactive material in front of 
the CAL and to particles lost in the rear beampipe, Wjb systematically underestimates the 
true W by approximately 13%, an effect which is adequately reproduced in the Monte Carlo 
simulation of the detector. Monte Carlo studies show a resolution of 6% for Wjb in the range 
120 < Wjb < 251 GeV. The reconstructed value Wjb is corrected to the true W by means of 
the Monte Carlo samples of events. 

The variables of the dijet system x® BS and [3 OBS are reconstructed as 

*SS = i ^-ffi'^ , and (3) 



JB 

Erpjet rf et 
jets ^T,cal e cal /.s 

Pcal 2x% l E n ' ' 



where x c is determined from the energies and angles measured in the CAL using 

/?2 _ p2 _ p2 _ p2 
v cal _ & r X r Y r Z (r\ 



7 



where Px (Py) is the X(y)-component of the vector P. Note that in the formulae for x® BS 
and (3 OBS many systematic uncertainties in the measurement of energy by the CAL cancel 
out. There are no significant systematic shifts in the variables x^f a f and f3° BS with respect 
to x° BS and [3 OBS . Therefore, no corrections are needed for these variables (x° BS ~ x^ B ^ 
and (3 OBS » /3°f 5 ). The resolution in x? BS ((3 OBS ) is 0.06 units in the region x° BS > 0.2 
(p OBS > 0.4). 

In the analysis presented here, large-rapidity-gap events are selected by using the vari- 
able rjmlxi which is defined as the pseudorapidity of the most-forward CAL cluster exceeding 
400 MeV The t] max variable for the Monte Carlo events at the CAL level is defined in the 
same way as in the data (i]™ l ax ). At the hadron level, t]^ x is calculated as the pseudorapidity of 



the most-forward particle with energy in excess of 400 MeV and pseudorapidity below 4.5 |17 
From Monte Carlo studies the resolution on rf^ ax is 0.1 units for v(^ ax < 1.8. 



7 Data selection 

Events from quasi- real photon-proton collisions were selected offline using criteria similar to 
those reported previously JI7J and briefly discussed here. A search for jet structure using the 
CAL cells is performed and events with at least two jets of E^* > 6 GeV and — 1.5 < rf et < 1 
are retained. The contamination from beam-gas interactions, cosmic showers and beam-halo 
muons is negligible after imposing a cut on the vertex reconstructed from three or more tracks. 
Neutral current DIS events are removed from the sample by identifying the scattered positron 



candidate using the pattern of energy distribution in the CAL |24 



Large-rapidity-gap events are selected by using the variable r(^ ax and events with r]™ ax < 1. 



are kept for further analysis [JPT| . The selected sample consists of events from e + p interactions 



with Q 2 < 4 GeV 2 and a median of Q 2 ~ 10 -3 GeV 2 . The event sample is then restricted to 
the kinematic range 134 < W < 277 GeV using the corrected value of Wjb- The data sample 
thus obtained consists of 403 events and represents 1.3% of the whole dijet sample in the same 
kinematic region except for the requirement on rf^ ax . The range in xjp spanned by the data is 
from 0.001 to 0.03 with a median of xjp « 0.009. 



8 Dijet cross sections 

Using the selected data sample of dijet events, differential dijet cross sections have been mea- 
sured in the above kinematic region with the most-forward-going hadron at r] max < 1.8. The 
cross sections have been measured as a function of r/- 76 *, E^ 1 , W, x° BS and f3 OBS for dijet 
production with E^ 1 > 6 GeV and —1.5 < rf et < 1. For each cross section, an integration over 
the remaining variables is implied. The cross sections refer to jets at the hadron level with a 
cone radius of 1 unit in the T] — <p plane. 

The Monte Carlo samples of events generated using POMPYT were used to compute ac- 
ceptance corrections to the r/ 6 ', E^, W, x® BS and (3 OBS distributions. These corrections take 
into account the efficiency of the trigger, the selection criteria and, the purity and efficiency of 
the reconstruction of jets. They also correct for the migrations in the variable f]^ l ax and yield 
cross sections for the true rapidity gap determined by rft^ x . The cross sections are obtained by 
applying bin-by-bin corrections to the distributions of the data. 
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8.1 Background and systematic uncertainties of the measurements 

The contribution from non- diffract ive processes has been estimated using a sample of direct 
and resolved processes generated with the PYTHIA Monte CarloQ. The fraction of large- 
rapidity-gap events in PYTHIA is strongly suppressed, although fluctuations in the final-state 
system may give rise to a rapidity gap in the forward region in a small fraction of events. The 
contribution from non- diffract ive processes as modelled by PYTHIA is approximately 20% in 
da/dW. This contribution, which has been subtracted bin-by-bin from the data, is listed in 
Table fj. 

Events with diffractively dissociated protons with masses less than ~ 4 GeV also con- 
tribute to the data set [|T^]. Since the pomeron flux factor of Donnachie and Landshoff accounts 
only for processes in which the proton remains intact after the collision and as the measure- 
ments are based upon a large-rapidity-gap requirement, the contribution to the measured cross 
sections from double dissociation has to be taken into account when comparing to model pre- 
dictions. A comparison of the measurements of the diffractive structure function in DIS based 
on the detection of the final-state scattered proton and on the Mjf-method [52], shows 



that double dissociation with M N < 4 GeV contributes (31 ± 13)% to the cross section. This 
estimation has been cross-checked with a study of exclusive vector-meson production |I3| yield- 
ing consistent results. Assuming Regge factorisation, the same contribution is expected to be 
present in the measurements of diffractive dijet photoproduction. Instead of subtracting the 
measurements for this contribution, the predictions of the MC program POMPYT (see next 
section) have been scaled up by the appropriate factor. Since this contribution affects equally 
the measurements of dijet cross sections and those of the diffractive structure function, the 
results of the QCD analysis (except for the overall normalisation) do not depend on its exact 
value. 

The statistical errors are indicated as the inner error bars in all the figures presented below. 
A detailed study of the sources contributing to the systematic uncertainties of the measurements 
was carried out [44]. These uncertainties are classified into six groups: 



The ?7mcLr variable in the data and simulated events was recomputed after removing the 
CAL cells with 7] > 3.25 in order to check the dependence on the detailed simulation of 
the forward region of the detector. This resulted in changes within ±10%. 

The energy threshold in the computation of ?7^L f° r data and simulated events was varied 
by ±100 MeV, yielding changes within ±10%. 

The amount of the non-diffractive contribution subtracted from the data was varied by 
±30%, yielding changes within ±15%. 

The relative contribution of quarks and gluons in the pomeron used in the simulated 
events was varied within the range obtained by the QCD analysis presented in the next 
section. The resulting changes are below 15% except for da/dx® BS in the lowest measured 
x OBS_ point (_|_4Q%) and da/df3 OBS in the highest measured (3 OBS -point (-20%). 

Variations in the simulation of the trigger and a variation of the cuts used to select 
the data within the ranges allowed by the comparison between data and Monte Carlo 
simulations yielded negligible changes. 



6 These calculations give a good description of the inclusive jet differential cross sections (without the large- 
rapidity-gap requirement) in the range — 1 < rf et < 1 flofl. 
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All these systematic uncertainties have been added in quadrature to the statistical errors and 
are shown as the outer error bars in the figures. 



The absolute energy scale of the cal jets in the simulated events was varied by ±5% H5j, 



resulting in changes of approximately ±25%. This uncertainty is the largest source of 
systematic error and is highly correlated between measurements at different points. It is 
shown SIS db shaded band in each figure. 

In addition, there is an overall normalisation uncertainty of 1.5% from the luminosity determi- 
nation which is not included. The results are presented in Table |l| 



8.2 Results 

The cross section da/drj^ et , shown in Figure 0, for the selected region of phase space has been 
measured in the rf et range between —1.5 and 1 integrated over Ej? > 6 GeV. The measured 
cross section dajdrf et shows a broad central maximum in if et . The cross section da/dE^ 1 
measured in the range of Ej 1 between 6 and 14 GeV and integrated over —1.5 < rf et < 1 is 
presented in Figure [5| The cross section do jdE^ 1 shows a steep fall-off as a function of E^ 1 . 
Note that in the data of Figures |2] and each of the two jets with highest E^ 1 in an event 
contributes to the cross section. 

The cross section da/dW for the selected region of phase space has been measured in the W 
range between 134 GeV and 277 GeV. The cross section, shown in Figure f|, falls for low values 
of W but remains fairly constant for large values of W . 

The cross section da/dx® BS has been measured in the x^ BS range between 0.2 and 1. The 
cross section da / dx° B , shown in Figure |5], peaks at high values of x° BS with a pronounced tail 
that extends to low-x° BS values. This result shows the presence of both a resolved- (low-x^ B ) 
and a direct-photon (high-rc^' 85 ) component in diffractive dijet photoproduction. 

The cross section da/d/3 , measured in the (3 OBS range between 0.4 and 1, is shown in 
Figure |6|. The cross section increases as [3 OBS increases and shows that there is a sizeable 
contribution to dijet production from those events in which a large fraction of the pomeron 
momentum participates in the hard scattering. 

The results are compared with the predictions of the MC program POMPYT using various 
parton distributions which are described in the following section. 



9 QCD analysis 

A combined QCD analysis of the ZEUS measurements of the diffractive structure function 
E®(j3,Q 2 ) in DIS [|13|] and of the measured dijet cross sections in diffractive photoproduction 
presented in the previous section has been performed following the proposal by Collins et al. |46| . 
This procedure assumes both hard-scattering factorisation and Regge factorisation of pomeron 
exchange in diffractive DIS and diffractive dijet photoproduction. 



The diffractive structure function F®(/3, Q 2 ) was measured [13|] for Q 2 values between 10 and 
63 GeV 2 , and (3 values between 0.175 and 0.65, by means of an integration over the measured 
range of xp, 6.3 • 10~ 4 < xp < 10 -2 . The fits to the DIS data are based on full NLO QCD 
calculations. On the other hand, the fits to the photoproduction data use calculations with LO 
matrix elements plus parton-shower as incorporated in POMPYT. Both direct and resolved 
processes are included. 
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In the calculations of the diffractive dijet cross sections using POMPYT, a s (^ 2 ) and the 
parton densities in the pomeron and the photon are evaluated at fi 2 = ft\. These computations 
may be affected by higher-order QCD corrections, which are expected to mainly change the 
normalisation (i.e. generate a X-factor). The agreement found between the PYTHIA calcula- 



tions of the inclusive jet differential cross sections and the ZEUS measurements |4(J indicates 
that in the case of the non-diffractive contribution the i^-factor is close to unity, within an 
uncertainty of ±30%. The i^-factor in the case of POMPYT is expected to be similar (with 
a similar uncertainty), as the same hard subprocesses are involved in the calculation of the jet 
cross sections. 

Each of the fits is represented by a parametrisation of the initial distributions at /i 2 , = 4 GeV 2 
for the quarks (f u /p(P, = fu/p = fd/p = fd/p) and for the gluon (f g /r). The other quark 
distributions are assumed to be zero at the initial scale. The parton distributions are evolved 
in /i 2 according to the DGLAP equations at NLO with the number of flavours set equal to five, 
the MS'-scheme with = 152 MeV and using the program from the CTEQ group ||47|| . 

The fits were performed simultaneously to the 1993 ZEUS measurements of F^(f3, Q 2 ) (Fig- 
ure |7|) and the measured cross sections da jdrf ei (Figure ^) and da/df3 OBS (Figure ||) in diffrac- 
tive photoproduction. These cross sections in diffractive photoproduction are those which are 
most sensitive to the shape of the pomeron parton densities. The following functional forms 
for the momentum-weighted parton densities in the pomeron at /x 2 , have been used in the fits: 

• hard quark + hard gluon: 

Pf u /p(P^l) = a 1 P(l-P) , Pf g/1P (P, t i 2 ) = b 1 p(l-P); 

• hard quark + leading gluon: 

PU /P (p,i4) = 0*13(1-13) , (3f 9/IP ((3,fi 2 ) = b 2 p 8 (l-(3) - 3 , 

• hard quark + (hard & leading) gluon: 

Pf u /p(P, l4) = 03/9(1 - (3) , PfMP, A*o) = h(3{l + c 3 p 8 (l - (3f 3 , 

where the values of Oj, 6, and Cj are determined from the fits. Note that all these functional 
forms are constrained to be zero at (3 = 1. The exact shape of the parton distributions are only 
weakly constrained for (3 > 0.8. 

The fraction of the pomeron momentum carried by partons is defined as 

Z% LO = f 1 df3 f3 [f g/F (f3,v 2 ) +Y,fu/A( 3 ^% 

JO j 

where the sum runs over all quark flavours which contribute at the given value of /i 2 . Since 
the pomeron is not a particle it is unclear whether or not the momentum sum rule (S^ LO = 1) 
should be satisfied. Therefore, S^ LO has been left unconstrained in the fits. 

The functional form of the "leading gluon" distribution is similar to one of the parametrisa- 
tions extracted by the HI Collaboration from the observed scaling violations of the diffractive 
structure function [|15|| . Parametrisations of the parton densities in which the pomeron is as- 



sumed to be made exclusively of quarks have been disfavoured by previous measurements [15, 17 



and are not considered here. Likewise, a parametrisation of the gluon density with a soft spec- 



trum has been ruled out [15, 17, 18 
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The fitted values of the parameters are shown in Table |^. These values have been obtained 
by subtracting from the data an estimated (31 ± 13)% contribution due to double dissociation. 
In the figures, the results of the fits have been scaled up to take into account this contribution. 
The results of the three different fits are compared to the measurements used for the fits in 
Figures || and [7[ The calculations based on these fits, which include a substantial hard 
momentum component of gluons at /z 2 , in the pomeron, give a reasonable description of the 
shape and normalisation of the measurements. The measured da /dj3 OBS rises as (3 OB increases 
while the calculations fall at f3 OBS > 0.8 for the parton distributions chosen. This comparison 
indicates that there is a sizeable contribution of gluons at large j3. 

The results of the fits have also been compared to those dijet cross sections in diffractive 
photoproduction which were not used in the fit (Figures § to y). The calculations provide 
a good description of the measured da/dE^ 1 , da/dW and da/dx° BS . The sensitivity of the 
measured da/dW to the value of a(0) in the pomeron trajectory is limited due to the cut on 
Eif 1 used to define the jets and the requirement on r] max . These comparisons are consistent with 
the conclusion that there is a large component of gluons with a hard momentum spectrum in 
the pomeron. Moreover, the measured da / dE^ 1 is well described by the calculations, indicating 
that the dynamics of dijet diffractive photoproduction is governed by the matrix elements of 
perturbative QCD. 

The predicted contributions from direct and resolved processes to the measured da/dx° BS , 
together with their sum, are shown in Figure || The sum of the contributions from resolved and 
direct processes gives a good description of the data. The shape of the contribution of either a 
purely direct-photon or a purely resolved-photon is not able to reproduce the data. A resolved- 
photon component is therefore needed in order to explain the shape of the measured cross 
section for values of x° BS below 0.8. This observation represents the first clear experimental 
evidence for the presence of both a resolved- and a direct-photon component in diffractive dijet 
photoproduction. Hard-scattering factorisation-breaking effects due to the resolved-photon 
component might be suppressed by the observed dominance of the direct-photon contribution 
to diffractive dijet photoproduction in the selected region of phase space. 

The fraction of the pomeron momentum carried by partons which is due to gluons, 

Zj jp J 

depends upon the scale at which the parton content of the pomeron is probed and has been 
computed for each of the fits. The determination of c^ LO (n 2 ) is affected by the following 
uncertainties: 

• The statistical and systematic uncertainties of the measurements used in the fits, which 
are the dominant sources of uncertainty. 

• The uncertainty (±30%) on the POMPYT calculations due to higher-order QCD correc- 
tions. 

• The uncertainty on the pomeron trajectory in the DL pomeron flux factor. The effect on 
c^ LO (fi 2 ) has been estimated by changing a(0) from 1.085 to 1.15. 

The central value of c^ LO (4 GeV 2 ) varies between 0.75 and 0.90 depending on the parametri- 
sation and is given in Table [|. Taking into account all the uncertainties of the three fits, the 
resulting range for c N g LO is 0.64 < c^ LO (4 GeV 2 ) < 0.94. 
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The fraction Cg LO {jj?) has also been computed at the /x 2 values probedf] by other measure- 
ments and is shown in Figure || a) /i 2 = (E 3 T et ) 2 = (6 GeV) 2 = 36 GeV 2 in dijet diffractive 
photoproduction (this analysis); b) /i 2 = (E^) 2 = (8 GeV) 2 = 64 GeV 2 in inclusive jet 
diffractive photoproduction 0; c) [i 2 = {E{ ) 2 = (20 GeV) 2 = 400 GeV 2 in dijet diffractive 
production in pp ]nj; and d) fi 2 = = 6400 GeV 2 in diffractive J¥-boson production in 
pp rF ' 



The result for c N g LO is consistent with the previous determination by ZEUS |T7] though now 
with an improved method and accuracy. The extrapolation of this result to the fi 2 values 
probed by the measurements in pp collisions is also consistent with the estimate of (70 ± 20)% 
given by the CDF Collaboration fll9|. However, care must be taken in these comparisons since 



the result presented here has been made using pomeron parton densities which depend upon 
the scale according to the DGLAP equations while those in |T7] and [T9| neglected any scale 
dependence. The use of different procedures does not have a significant effect on the extracted 
values of c g , but affects substantially the values of Ejp. 

We conclude that for the selected region of phase space, which includes the r\ max requirement, 
and within our experimental uncertainties it is possible to reproduce the measurements of the 
diffractive structure function JL3| and of the dijet cross sections in diffractive photoproduction 
with the Ingelman-Schlein model PJ. In this model, which assumes Regge factorisation, the 
parton densities of the pomeron evolve according to the DGLAP equations. The dominance of 
the direct process for the measured region of phase space could be limiting the observation of 
hard-scattering factorisation-breaking effects due to the resolved-photon component. The data 
require the fraction of the pomeron momentum carried by partons which is due to gluons to lie 
in the range 0.64 < cf LO (4 GeV 2 ) < 0.94. 



10 Summary and conclusions 

Measurements of the differential cross sections for dijet photoproduction with a large rapidity 
gap in e + p collisions at a centre-of-mass energy of 300 GeV have been presented. The e + p 
dijet cross sections refer to jets at the hadron level with a cone radius of one unit in the 

11 — ip plane. They are given in the kinematic region defined by Q 2 < 4 GeV 2 (with a median 
Q 2 « 10 -3 GeV 2 ) and 134 < W < 277 GeV with the most-forward-going hadron at r] max < 1.8. 
These cross sections have been measured as a function of rf et , E 3 ^ 1 and W for dijet production 
with E J T et > 6 GeV and -1.5 < rf et < 1. 

The measured cross section da/dx° BS as a function of x° BS , the fraction of the photon 
momentum participating in the production of the two jets with highest Ej?, peaks at x° BS ~ 1 
with a pronounced tail to lower values. This result is clear evidence for resolved- and direct- 
photon components in diffractive dijet photoproduction. 

A measurement of the cross section for diffractive dijet photoproduction as a function of 
j3 OBS , the fraction of the pomeron momentum participating in the production of the two jets 
with highest E^ 1 , has been presented. For the selected region of phase space the measured 
cross section da /d(3 OBS increases as (3 OBS increases. This result shows that there is a sizeable 
contribution to dijet production from those events in which a large fraction of the pomeron 
momentum participates in the hard scattering. 



A QCD analysis of the measurements of the diffractive structure function in DIS [T3[ and 



7 For the non-DIS measurements referred to here the /i 2 value is not uniquely defined and has been 

computed at either (E^) 2 or the square of the W-boson mass depending on the specific final state. 
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of the measured cross sections dajdrf et and da/df3 OBS presented here has been performed. 
The pomeron is assumed to have hadron-like partonic structure in the form of parton densities 
which evolve according to the DGLAP equations. It is possible to reproduce both sets of 
measurements when a substantial hard momentum component of gluons in the pomeron at the 
initial scale of 2 GeV is included. The data require the fraction of the pomeron momentum 
carried by partons which is due to gluons to lie in the range 0.64 < c^ LO (A GeV 2 ) < 0.94. 
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Table 1: Measured cross sections in the kinematic region defined by Q 2 < 4 GeV 2 and 
134 < W < 277 GeV with the most-forward-going hadron at r] max < 1.8. The contribution 
from non-diffractive processes, which is given in the last column, has been subtracted. The 
measurements contain an estimated (31 ± 13)% contribution from double dissociation. The 
statistical and systematic uncertainties —not associated with the absolute energy scale of the 
jets— are also indicated. The systematic uncertainties associated to the absolute energy scale 
of the jets are quoted separately. The overall normalisation uncertainty of 1.5% from the 
luminosity determination is not included. 
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Table 2: Fitted values of the parameters for each of the three fits discussed in the text. The 
values of c^ LO (A GeV 2 ), S^ LO and xltati an d the number of degrees of freedom (dof) for each 
of the fits are also shown. 
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Figure 1: Examples of Feynman diagrams for diffractive dijet photoproduction: a) direct 
process and b) resolved process. 
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Figure 2: Measured jet cross section da ldrf et in dijet events (see text) integrated over > 
6 GeV in the kinematic region defined by Q 2 < 4 GeV 2 and 134 < W < 277 GeV with the 
most-forward-going hadron at r] max < 1.8 (black dots). The contribution from non-diffractive 
processes (see Table [L|) has been subtracted. The measurements contain an estimated (31±13)% 
contribution from double dissociation. The inner error bars represent the statistical errors of the 
data, and the outer error bars show the statistical and systematic uncertainties —not associated 
with the absolute energy scale of the jets— added in quadrature. The shaded band displays 
the uncertainty due to the absolute energy scale of the jets. For comparison, the results of the 
QCD fits, which have been scaled up to account for the contribution from double dissociation, 
are shown (see text). The results of the QCD fits have been obtained by an integration over 
the same bins as for the data and are presented as smooth curves joining the calculated points. 
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Figure 3: Measured jet cross section da/dE^ 1 in dijet events (see text) integrated over —1.5 < 
rf 6% < 1. Other details as in Figure |2|. 
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Figure 4: Measured dijet cross section da/dW integrated over E 3 ^ > 6 GeV and -1.5 < 
rjjet < i Other details as in Figure |2|. 
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Figure 5: Measured dijet cross section da/dx® BS integrated over E 3 ? > 6 GeV and -1.5 < 
rjjet < 1. Other details as in Figure |2|. For comparison, the calculations for the resolved (dot- 
dashed line), direct (dashed line) and resolved plus direct processes (solid line) based on the 
QCD fit with the "hard quark + hard gluon" parametrisation (see text) are shown. 
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Figure 6: Measured dijet cross section da/df3 OBS integrated over E 3 T et > 6 GeV and -1.5 < 
rjjet < i Other details as in Figure 0. 
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Figure 7: Measurements of F®(j3,Q 2 ) in DIS @ as a function of (5 for fixed values of Q 
compared to the results of the QCD fits: "hard quark + hard gluon" (solid lines), "hard 
quark + leading gluon" (dashed lines) and "hard quark + (hard & leading) gluon" (dotted 
lines) parametrisations. The measurements contain an estimated (31 ± 13)% contribution from 
double dissociation and the results of the QCD fits have been scaled up to take into account 
this contribution. 
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Figure 8: The values of c^ i0 (/i 2 ) as a function of /i 2 for each of the QCD fits described in the 
text. For the "hard quark + hard gluon" (upper band) and "hard quark + leading gluon" (lower 
band) parametrisations, the central values (indicated by the lines), the statistical and systematic 
uncertainties —not associated with the absolute energy scale of the jets— added in quadrature 
(light-shaded bands) and the uncertainty due to the absolute energy scale of the jets (dark- 
shaded bands) are shown. The theoretical uncertainties are included in the light-shaded bands 
and discussed in the text. For the "hard quark + (hard Sz leading) gluon" parametrisation, 
only the central values are shown (dotted line). The /i 2 values at which (/z 2 ) has been 
computed are indicated by the arrows. 
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